Universal cumulants of the current 
in diffusive systems on a ring 
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We calculate exactly the first cumulants of the integrated current and of the activity (which is the total number 
of changes of configurations) of the symmetric simple exclusion process (SSEP) on a ring with periodic bound- 
ary conditions. Our results indicate that for large system sizes the large deviation functions of the current and 
of the activity take a universal scaling form, with the same scaling function for both quantities. This scaling 
function can be understood either by an analysis of Bethe ansatz equations or in terms of a theory based on 
fluctuating hydrodynamics or on the macroscopic fluctuation theory of Bertini, De Sole, Gabrielli, Jona-Lasinio 
and Landim. 
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I. INTRODUCTION 

The symmetric simple exclusion process (SSEP) | 
is one of the simplest lattice gas models studied in the theory 
of non-equilibrium systems. It consists of hard-core particles 
hopping with equal rates to either of their nearest neighbor 
sites, on a regular lattice. At equilibrium, when isolated, the 
system reaches in the long time limit an equilibrium where all 
accessible configurations are equally likely. Also, when equi- 
librium is achieved by contact with one or several reservoirs 
at a single density p, all sites are occupied with this density p 
and the occupation numbers of different sites are uncorrelated. 

As soon as the system is maintained out of equilibrium, by 
contact with reservoirs at unequal densities, there is a current 
of particles and one observes long range congelations in the 
steady state fs*]. In this out of equilibrium case several ap- 
proaches have been developed to calculate steady state prop- 
erties, such as the fluctuations or the larg e deviations of the 
density or of the current H S i H [H [H d [H [H [B [H 

Wim- 

A lot of progress has been made over the last decades on 
the study of the fluctuations and the large deviation func- 
tions of the current in equilibrium or non equilibrium sys- 
tems. The large deviation function of the cuiTent can be 
viewed as the dynamical analog of a free energy, as dis- 
cussed by Ruelle in the early seventies ifToll . The idea back 
then was to build up a thermodynamic formalism based upon 
probabilities over time realizations rather than over instanta- 
neous configurations. Generic properties of these large de- 
viation functions were later discovered such as the fluctua- 
tion theorem which determines how the large deviation func- 
tion of the current is changed under time reversal symmetry 

s III 111 m m m Hi [23, m . 

In the present work, we obtain exact expressions for the first 
cumulants of the integrated current and of the activity (which 
is the number of changes of configurations) during a long time 



t for the SSEP consisting of N particles on a ring of L sites. 
For large system sizes, these cumulants and the associated 
large deviation functions take universal scaling forms. We 
show how these scaling forms can be calculated for the SSEP 
by the Bethe ansatz or for more general diffusive systems on 
a ring by a theory based on fluctuating hydrodynamics or on 
the macroscopic fluctuation theory developed by Bertini, De 
Sole, Gabrielli, Jona-Lasinio and Landim M |Tol Hej ITtI ITsIl . 
In the Bethe ansatz approach these scaling forms can be ex- 
tracted from a detailed analysis of finite size effects similar to 
what was developed recently for quantum spin chains in the 
context of string theory ll29[|30ll . In the fluctuating hydrody- 
namics approach, it results from the discreteness of the wave 
vectors of the fluctuating modes on the ring. 

Universal distributions of the cuiTent characteristic of the 
universality class of the KPZ (Kardar-Parisi-Zhang) equation 
jfll m, |33l H, have been calculated in the past I35l l36i IstI 
l38ll for the asymmetric exclusion process (ASEP). The distri- 
butions obtained in the present paper are different and belong 
to the Edwards-Wilkinson universality class ll39ll . 

We begin by presenting in Sec|ll] exact expressions of the 
first cumulants of the cuiTent and of the activity for the SSEP 
on a ring . This is where we see that the cumulants of the 
integrated cuiTent and of the activity take scaling forms when 
the size of the ring becomes large and where emerges the idea 
that the large deviation function of the current and of the ac- 
tivity obey the same universal scaling function. This is con- 
firmed in Sec|lII|by Bethe ansatz calculations. By resorting to 
fluctuating hydrodynamics in SeclIVIwe are able to formulate 
the particular case of the SSEP within a more general frame- 
work using the Bertini, De Sole, Gabrielli, Jona-Lasinio and 
Landim approach and to show that the same universal distri- 
bution of the current fluctuations of the current are present in 
a larger family of diffusive systems. 
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n. EXACT EXPRESSIONS OF THE FIRST CUMULANTS 

We consider a system of N particles on a one-dimensional 
lattice of L sites with periodic boundary conditions. Each site 
is either empty or occupied by a single particle. A microscopic 
configuration C = {ni}i,i....,L can be specified by L occupa- 
tion numbers (where ji^ = 1 if site i is occupied and = 
if site i is empty). In the simple symmetric exclusion process, 
SSEP, each particle hops to its right neighbor at rate 1 or to 
its left neighbor at rate 1, provided the target site is empty. 
In the present paper we try to determine the distribution of 
the total integrated current Q{t) and of the total number KU) 
of changes of configuration (that we will call the activity [4011) 
during a time interval (0, i). To do so we define the generating 
functions of the cumulants of Q and K as 



lim 

t — *C30 



t 



ln(e--^) 
TpKKs) = lim , (1) 



where the brackets denote an average over the time evolutions 
during the time interval (0,<). As the evolution is an irre- 
ducible Markov process with a finite number of states, the 
long time limits in ([T]i do not depend on the initial configu- 
ration and the generating functions defined in ^ can be cal- 
culated as the largest eigenvalue of a matrix 

Because the calculations are very similar for both observ- 
ables K and Q, we shall first focus on the activity K and 
explain how to calculate the cumulant generating function 
iPk{s) as a perturbation series in powers of s. We will then 
present only the results for iIjq{s). 



A. The cumulants of the activity K{t) 

In order to determine t/^x, as in [36], one can write a Mas- 
ter equation for the probability P{C, K, t) to find the system 
in configurations C at time t, given that the activity at time 
t is K (i.e. given that the system has changed K times of 
configurations during the time interval (0, t)). 

dtP{C, K, t)^-r{C)P{C, K, W{C' C)P{C',K^1, t) 

c 

(2) 

where W{C C) is the transition rate from configuration C 
to C, and r(C) W{C — > C) is the escape rate from 

configuration C. 

If one introduces the generating function P(C, s, t) — 
e^'*^P(C, t), its evolution satisfies 

dtPiC, s,t)=J2 C')P{C', s, t) (3) 



where 



IK 



{C,C')=t-'W{C' ^C)-r{C)5cfi' 



(4) 



In the long time limit, P{C, s, t) grows (or decays) expo- 
nentially with time, with a rate given by the eigenvalue with 
largest real part |36] of the modified matrix Wx- Thus iPk{s) 



can be calculated as this largest eigenvalue of . For s = 0, 
Wk reduces to the evolution operator of the Master equa- 
tion W for the symmetric simple exclusion process, and this 
largest eigenvalue (which is 0) as well as the related eigenvec- 
tor are known. We now present a way of obtaining the large 
deviation function tpK, by a perturbative expansion ||4TI 
in powers of s. 

The idea is to start from the eigenvalue equation for ipK and 
its eigenvector P, 



i^K{s)P{C, 



c 



k{C,C')P{C',s) 



(5) 



normalized such that P{C, s) — 1. One can then define 
the average {A{C))s of an observable .4(C) in the correspond- 
ing eigenstate, (i.e. {A{C))s = X^c -^(^)^(^' 
is the same as averaging, in the limit of a long time interval 
(0, t), over all trajectories weighted by a coefficient e^'*^'*'). 
Note that, though the value of K{t) is defined on trajectories 
running from to t, the observable A{C) is evaluated at the 
final time t. From the eigenvalue equation, one gets 

^k{s){A{C))s = e-^ W{C ^ C')A{C')^ -{A{C)r{C)), 

(6) 

where the escape rate r{C) is twice the number of clusters of 
adjacent particles in the system 



r(C) = ^ W(C ^ C) = 2^n,(l - n,+i) . (7) 

Choosing A{C) = \ in^ leads to 

^k{s) = (e^^ - l)(r(C)), = 2L(e-^ - l){p - C7,(l)) (8) 

where Cs{r) — {nini^r)s is the correlation function (which 
by translational invariance does not depend on i) computed 
within the eigenstate P(C, s), and p — N/L is the average 
density. 

For the leading contribution as s ^ 0, we can use the fact 
that at s = the eigenvector is known (this is the equilibrium 
distribution, for which all allowed microscopic configurations 

are equally Ukely), so that '4>k{s) = -'2.N {l - ^^fi^ s + 

0{s^). In order to compute the O(s^) contribution from ([8]l, 
we need to evaluate Cs(l) at order s, which can be done by 
choosing A{C) = riinj in (|6]l. This requires the knowledge 
of the correlation function Cs(r) = {nini+r)s at order 0{s). 



For .4(C) 



in (|6]l one gets 



Cs{1)-C,{2) = sAn,l + 0{s^) 

N{N - 1){L - N){L - N-1) 



where Ajyx 



L{L-lf{L-2) 



2A 



N.L 



+0{s^) (9) 



L-3 
for 2 < r < L - 2 
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which have the following solution 

N{N-1) 6r{L-r)-L{L + l) , , 



L(L-l) 



6(L-3) 



(10) 

We can therefore extract t/jK up to ©(s^) and {K'^)c/t fol- 
lows. 

To obtain higher cumulants, we have repeated the same pro- 
cedure, with the observables -4(C) = niUjUk and -4(C) = 
niUjUkni. The calculations are longer but very similar We 



found that the first cumulants of K, lim(^oo(^")c/^ 



(-1)" tI^ 

L and of 



s=0 



, when expressed in terms of the system size 



a{p) = 2p(l - p) 



are given by (in the t oo limit) 



2N{L-N) 
I? 



(11) 



L^a{L^a + 4L - 4) 



t 6(L-1)2 
{K^)c _L^(t[- L^<7^ + L^a{2 + 3a) - 2L^g + 48(L - 1)^] 



t 



60(L - If 



= L^(y[a*L''{lQL'' - 70L^ + 175L - 153) - 4cr^L'*(L - l){nL^ - 69L^ + 154L - 126) 
^16aL^{L - 1)2(3^3 - 17L^ + 46L - 63) + 2112(L - lf{L - 3)) (2520(L - 1)*{L - 3)) \ (12) 



When L becomes large, while p — N/L is kept fixed, the 
asymptotic behavior of the above cumulants reads 



t 



t 



aL, 
60 ' 



t 



(13) 



t 



252 



One might have expected the derivatives at s = of the eigen- 
value ipK to become extensive for a large system size L (after 
all, as we shall see it in section |III] it is always possible to 
view as the ground state energy of a short range Hamil- 
tonian). Yet this is not the case since the second and higher 
cumulants grow faster than linearly with L at fixed density 
p. This suggests that, in the large L limit, i/jk / L becomes a 
singular function of s at s = 0. 

Also one can guess from ST3[ that for n> 2 

nr2n~2 

^ a L 

t 

oo and s ^ 0, the eigenvalue ipK takes a 



and that for L 
scaling form 

lim 



(14) 



where the scaling function Tk is given by 

Tk{u) = \u^ + ^u^ + ^u* + 0{u''). (15) 

We shall see in sections |III] and |IV] that this scahng function 
can be fully determined and written as 



(u) — —4 ^nn \/n^TT^~-2' 



u — t^^-n^ + u 



(16) 



or equivalently (see appendix A) as 



B 



2k-2 



k>2 



(fc- 1)! k\ 



{-2uY 



(17) 



where the Bernoulli numbers i?„ are known to be simply the 
coefficients of the expansion x{e^ — 1)^^ = BnX^ /n\. 
As a consequence, the generalization of ( fT3] l will be for n>2 



B2n-2 nr 2ri-2 
-a L 



{n-l)\ 



(18) 



B. The cumulants of the current 



The same procedure can be followed for the total integrated 
current Q (which can be defined by Q = X]j=i where 
Xj (t) is the total displacement of the jth particle during the 
time interval (0, t)). Its cumulant generating function ipQ de- 
fined in ([T]i is the eigenvalue (with largest real part) of the 
matrix 



Wq(C,C') = W{C' C)s-''^^''^'> - r{C)dc,c 



(19) 



where j{C' ,C) is +1 or —1 depending on whether a particle 
has moved to the right or to the left when the system jumps 
from configuration C to configuration C. Using an expansion 
in powers of s as in III A I we have obtained (in the limit t ^ oo) 



4 



L^a'^ ((10L4 - 2^3 + 27L2 - I5i + is)^^ - 4(L 



1 



4^2 



L-1' < 2(i-l)2 

LV2 ((L2 - L + 2)cr - 2(L - 1)) 



4(L- 

l)(lli2_L^ 



-l)3(L-2) 
12)cr + 48(i 



1) 



24(L- - 2)(i - 3) 



(20) 



(21) 



with the corresponding large L behaviors (for p — N/ L fixed) down otherwise. In this basis one finds that 



t 



t 



4 {Q')c 



- —7-2 

- 2^ ' 



t 4 ' t 12 

As for K, these results indicate that for n > 2 



(22) 



n 7-2n-2 



and that tjjq takes a scaling form, in the limit L 

s ^ 



lim 

L — >oo 



.2 (Q2 



-2.252 



00 and 



(23) 



where, according to (|22]) . the expansion of JTq (u) in powers 
of u coincides with the expansion ( fTSl l of !Fk{u), at least up 
to the 4th order in u. 

We will see, in section HVl that these two scaling functions 
(which appear in (fT4l) and in ( |23] )) are in fact the same. There- 
fore the formula which generalizes ( |22] | will be for n > 2 



(2n)! B2„_2 
2" (n- 1)! n! 



r 2n— 2 
(7 L 



(24) 



III. BETHE ANSATZ 

It is well known that the Bethe ansatz allows one to cal- 
culate the eigenvalues of matrices such as Wk(C,C') and 
WqiCX!) defined in dgjlgll for exclusion processes il M, 
[33MlillillillillM|4iR ^,'5^. I n this section we show 
how to obtain the scahng forms (I14I23I) from the Bethe ansatz 
equations. 



-Wk 



^Q = ^-\il [cosh, + afa,^i) + afa^+i(25) 



1=1 

-zsinhs«fT^i-afa^i)] ===-Wq 

where we have resorted to the Pauli matrices cr^'^'^. In this 
language, the quantities ipx and t/jq are the ground state en- 
ergies of these operators. It also suggests that the methods of 
one-dimensional exactly solvable models apply in our case, 
such as the Bethe ansatz, as was exploited for similar systems 
in the past L36„.52il53|]. 

As the number of particles on the ring is fixed, we 
need to find the ground state with a fixed particle density 
p, that is, at fixed transverse magnetization J^i'^i- The 
quantum operators appearing in ( |25] ) have of course been 
extensively studied ll54ll . including within the framework 
of stochastic dynamics ifssll . For instance, following the 
notations of Baxter ifsill the operator e^Hx: is the ferro- 
magnetic XX Z chain with anisotropy parameter A = e*. 
Similarly, Hq corresponds to an XX Z chain with additional 
Dzyaloshinskii-Moriya interactions. A study of an operator 
closely related to Hq was carried out by Kim in 1995. 
His results will be recalled at the end of the present section. 

The Bethe ansatz consists in looking for the ground state 
of Hk or Q in the form of a linear combination of -particle 

plane waves (see (4^ I52I1 ). We denote by {xj}j^i n the 

positions of the TV particles and we postulate that the right 
eigenvector of Wk can be cast in the form 



N 



(26) 



A. Relation to spin chains 

It is possible to write the matrices WxiCC) and 
Wq{C,C') as quantum spin-chain Hamiltonians IsTf]. We use 
the correspondence in which the z component of a two state 
spin operator is up when a particle is present at site i, and is 



where V = (p(l), • • • ,p{N)) is a permutation over the first 
N integers, and the Q's are a priori complex numbers. This 
is an exact eigenstate provided these parameters satisfy the 
so-called Bethe equations. These take different forms for K 
and Q. We now discuss how to implement the Bethe ansatz to 
calculate tpK{s) and tpgis) defined in ([l). Technical details 
have been gathered in the appendices. 
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B. Bethe ansatz for K 



and 



For the expression (l26b to be an eigenvector of Hk or Wg 
the Cj 's have to satisfy a number of constraints |56], the so- 
called Bethe (see for example |49]) equations 



TV 

n 



1 - 2e^0 + CO 



1 - 2e^Cj + GO 



(27) 



The expression of iPk{s) is given by 



AT 



27V 



(28) 



Our goal is to obtain (fT4l) from dZTl i and (l28T l in the double 
limit s ^ and L ^ oo keeping sL^ and N/ L — p fixed. 
Because of the particle-hole symmetry the discussion below 
is limited to the case p < 5- 

In the large L limit, the (j 's accumulate on a curve which 
depends on s and as s — > 0^ becomes a finite arc of the unit 
circle (see [^,'56'] and references therein). Note however that 
the s > case can be approached by similar methods. 

If one writes 



cos 5 



(29) 



becomes 



1 ^ 



where 



t/(fci, kj) ^ — In 



1 - 2e*'='^ cos (5 + e*('=^ 



1 - 2e*'^3'5 cosJ + e*('=j+'=')'5 



(30) 



(31) 



(32) 



In the limit 6^0, one can check that when ki — kj — 0{1) 



U{h,k,)^2\-^+0{S') 



ki kj 



(33) 



In the large L limit, however, the distance between consecu- 
tive ki becomes of order 1/L ^ S and for i~ j of order 1 one 
should use instead 



U(^ki^ kj 

Therefore one can rewrite ( [3T| i as 



10 



ki - kj + iS{\ - kf) + i5ki{ki - kj) - 5'^ki{\ ~ kf) 



16(1 - ki) - iShih - kj) + 6^hil - kf) 



Lki 



-"o<j<*+no 



fcj - kj + i6{l - kf) + iSki{k, - kj) - 6^k^{l ~ kf) 
ki - kj - iS{l - kf) - i6h{ki - kj) + S^h{l - kf) 



+ 



E 

j^li-no,i+no] 



(34) 



(35) 



where no is a fixed large number 1 <C n-o ^ i, so that one 
can use expression ( |33] l for \ j — i\ > and ( [34l l for \i — i\ < 
riQ. As shown in appendix B, the two sums (1891 1001 ) in ( |35] | 
depend on the cut-off no but this dependence disappears when 
the two terms in the right hand side of (1351) are added. 

In the large L limit, the ki become dense on an interval 
(—0, 0) of the real axis, with some density g{k). In what fol- 
lows we will assume that the ki are regularly spaced according 
to this density, meaning that 



L j g{k)dk = j-i and L j g{k)dk = N . (36) 



Replacing the two sums in ( [35] ) by their expressions ( 18911 OOl l 
obtained in Appendix B, one gets that for k — ki the density 



g{k) should satisfy 



k = 2V I Ak'g{k') 
J-e 



k-k' L 



g'{k){i-k^) 
g{k) 



2k 



X7r(l - k'^)g{k)LS coth[TT{l - k^)g{k)L5] 

(37) 

If we make the change of variable k' = Oy, k = 9x, and 

g{k){l - e) = ^{x) (38) 
equation (l37l l becomes 



V / dy- 



1 y-x 



fix) 



(39) 



6 



where 



and 



f{x)=-— + ^ _Z^l^^coth[L^7r0(x)] + ... 

i— ,— . ^^^^ 
As explained in ( I101I102| | of Appendix C one can invert 

and express (f){x) in terms of f{x) 



(47) 



C 



7r2yr 



1 2^ 



where the constant C is so far an arbitrary constant. 
For small 6, one can write (l28T l. using ( I30I38I41T |. as 



(41) 



So far, the constant C remains undetermined. 

The leading order corresponds to using expression ( [33] ) in 
OTT l even when i and j differ by a few units. For the contin- 
uum description to be valid, we are now going to argue that 
(j>{x) should remain finite as x — > ±1, or, in terms of the orig- 
inal density g, that g{k) remains finite as A: ^ ±6*. This will 
impose (see ( |44] |) that 



C 



iPk{s) ~ V^^(l - kf) ~ LS^ / g{k){l - k^)dk 



in 



dx- 



C 



1 \/l — x"^ 7r2\/l — x'' 



1 y-x 



which gives using ( I113I115I ) 
Also, as (l36ll 



(42) 



g{k)dk = p 



one has (ISsBTT i 



p = 9 I dx 



C 



(l-02^2-)yf 



7r2(l- 6*22:2) Vr 



V i' ^^^^fiy)dy 



-1 y-x 



which can be simplified using ( II 1611201 ) 



P = 



dy. (43) 



C. The leading order in tlie large L limit 

For large L (at fixed LS), (|40| | reduces to f{x) - 
so that dm becomes to leading order using ( IllOl l 



47r\/l — a; 



whereas ( |43] l becomes using ( 11201 ) 

COir 1 6l2 - 2 



P = 



-6ia;/2, 
(44) 

(45) 



Therefore for a fixed density p of particles, the constant C 
in ( 14 11441 ) and the eigenvalue ( l42b are given, to leading order 
in -i, by 



C ^ — 



- 1] Vi-0^ + 



2- 



(46) 



Indeed if we order the N solutions ki and focus on the ones 
closest to 6*, . . . < k^^i < k^ < 9, then we may estimate 
using ( [36] l the difference between k^ and 9, or between fcjv_i 
and kN- If C ^ then g{k) ^ {9 - k)'^/^ as k ^ 9 
implies that fcjv — ^jv-i ~ L^^ . This is not compatible with 



kN > 



(which follows from (I31I33I) ). where the 



L A; TV — ^TV — 1 

right hand side of this inequality would be 0{L) in contradic- 
tion with the fact that k^ < 9. Hence we must have AttC = 9, 
in which case k^ ~ kjsi^i ^ i^2/3 ^jj^j (jiei-g is no contradic- 
tion. 

It then follows that 



and therefore tpK(s) 



' - Vp(i - p) 

L5^p{l~ p) anddH 



9VT 



0(i 

27r \L 



D. The next order 



(48) 



(49) 



Once (f) is known to leading order ( |49] i, one can update the 
expression (00]) 



9x (1-02^ 
fix) = d cotn 



LS9VT 



and one gets from (143] ) 



C9t: 



1 



47r\/l - 612 7_i 



y^coth ^ — dy 



(50) 



(51) 



Then using the fact that (see (l80] ) in appendix A) 



we get 



VT~9^ 2 4\/T~~e2 4LVT~~^ 



1 



^252^2 



(53) 



8 



7 



and this gives ( |42] | 

V'K (s) = L<52C7r6l = L^^ 



2 - 



02 



1 



4L ^ L3(52-^ ' 



4 
8 



The leading order (the first two terms of d54] i) has a mini- 
mum for 9 given by (|48] |. Therefore to obtain iPk{s) at first 
order in one can simply replace 6 by (l48T l in ( l54b and one 

gets 



1 



1 



4 V 

which is equivalent (see (129148b ) to (fT4l i 



1 / L2,5202 



(55) 



It is shown in ( [85] ) of appendix A that for large negative 



replaced by 



N 

n 



e^ - 2C. 



e« - 20 



^00 



(59) 



(54) Given the solutions Q to ( |59] l, the expression of t/iq reads 



-2Ar + e-"[Ci + ... + Cw] +e" 



1 1 
- + ... + — 

Ci 

(60) 

By a method following closely the steps of the Bethe ansatz 
for K, the basic ingredients of which are provided in appendix 
E, we arrive at the following result for tpq, 

i^Q{s) = laip)s'{L + l) + L-'T^-^^^^^ (61) 
which leads to the asymptotic behavior as L ^ oo, 



27/2 

Stt 



(56) 



This implies that ( fT4b becomes for small negative s (but large 
negative L^s) 



(57) 

So for s small, but L^s large, the extensivity of iPk{s) is re- 
covered and dSTl l gives the beginning of the small s expansion 
in the large L limit. 

One can also notice that the function J-{u) (fT6] l becomes 

2 

singular as u ^L. xhis indicates the occurrence of a phase 

transition discussed at the end of section |IV| for u > ^ the 
optimal profile to reduce K is no longer flat and the system 
adopts a deformed profile as in lfl6ll . In fact in the limit s 
+00 the configurations which dominate are those formed of a 
single cluster of particles and the activity is limited to the two 
boundaries of this cluster. 

The result ( fSSt or equivalently ( Il4] i with T given by ( fTSI l 



n>l 



2u- 



2 2 
n TT 



(58) 



gives the leading finite-size correction to ipxis). These finite 
corrections have been calculated recently, starting from the 
Bethe ansatz equations, for several spin chains in the context 
of string theory and expressions very similar to our T have 
been obtained [29.1. Note also that a more systematic approach 
has been developed to calculate the next finite size correction 



E. Bethe ansatz for Q 



L Z OTT 



(62) 



The Bethe equations ( |59] l are very close to that considered 
by Kim |52] who worked out the asymmetric exclusion pro- 
cess case. As outlined in appendix E, it seems that Kim's 
results cannot be extended to the SSEP. We think that this is 
at the origin of the discrepancy between our expression (l62b 
and what was found earlier (expression (A. 12) of ll20ll ) for the 
same quantity ijjQ{s). 

Before concluding this section devoted to the Bethe ansatz, 
let us mention that, both for the current or the activity, one can 
obtain ipq (s) or ipK (s) in the s oo limit by directly solving 
(|59] l or ( |27| i. We do not give these expressions here because 
they are out of the universal regime. 



IV. FLUCTUATING HYDRODYNAMICS AND THE 
MACROSCOPIC FLUCTUATION THEORY 

In this section we are going to show that the expressions 
(I14I23I I can be recovered by a macroscopic theory based on 
hydrodynamical large deviations |[ll0,01- 



A. Calculation of ipg for a general diffusive system and 
derivation of f23l ) 

The macroscopic fluctuation theory developed by Bertini, 
De Sole, Gabrielli, Jona-Lasinio and Landim |6, 7, 8, 9, 10] is 
based on the fact that, for a large system of size L, the density 
and the current of a diffusive system take scaling forms. If one 
defines pi{t), the density averaged in the neighborhood of site 
i at time t, and Qi{t), the total flux between site i and i + 1 
during time t, these quantities take scaling forms iflTlfTsIl 



The eigenvector corresponding to the largest eigenvalue of 
Wq can be written as in ( l26] l, with the Bethe equations (l27l i 



(63) 



8 



Q^it) = LQ 



t 



This allows one to define a rescaled current j{x, r) as 



(64) 



dQ{x,T) 



The average microscopic current between site i and i + 1 is 
related to the rescaled current j by 



dt 



1 



t 



From the macroscopic fluctuation theory 10, 0. H. [13. fisll . 
the probability of observing a rescaled current t) and a 
density profile p{x, r) over a time t ~ TL^ is given by 



exp 



.T .1 

-LI dr dx 
Jo Jo 



Pro({p(x,T),j(a;,r)}) ' 

[jix,T)+Dipix,T))p'ix,T)]^ 

2a(p(x,T)) 



(65) 



where the current j{x, t) and the density profile r) sat- 
isfy the conservation law 



dp 
d^ 



dx 



(66) 



and the diffusive system under study is characterized by the 
two functions D{p) and (j{p). For the SSEP, these functions 
are known: D{p) = 1 and a{p) = 2p{l — p) (see j^])- 

Note that ( |65] | can be seen as the fact that the macroscopic 
density p{x, r) and the macroscopic current j{x, r) satisfy in 
addition to the conservation law (|66T l a Langevin equation of 
the form lUt]. 

j{x, t) = -dxp{x, t) + £,{x, t) (67) 
where f (x, t) is a Gaussian white noise 

{ax,T)i{x\T')) ^ L-^a{p{x,T))6{x-x')5{T~T') . (68) 

The contribution of a small time dependent perturbation to 
a constant profile po and a constant rescaled current jo^ 

p(a;,T) = po + Sp{x,T) 

j{x,T) = jo + Sj{x,T) 

to the quadratic form in ( |65] l is 

[j(x,t)+£>(p(rr,t))p'(x,t)]^ ^ Jo_ Jo . . _ 

2(7(p(a;,t)) 2a a ^ 2a^ ^ 

joD , 5f + 2DSjdp' + i^^^p'^ + 2joD'Sp6p' 

+ dp + 

a 2(7 



joa'{Sj6p + DSpSp') 



(69) 



where the functions D,(j,a, a" are evaluated at the density 

Po- 

If one considers a fluctuation of the form 



— ibjt—ihx 



(70) 



one has 



6p' = i/c [afe^j^e 



ai. „,e 



and due to 



Sj = -cj[ 



/ il^ iujT+ikx I * — i(i;t— ifcx] 



The ring geometry (x = x+1) imposes that the wave numbers 
k are discrete 

k = 27rn with n > 1 

Also because one considers a finite time interval T, the fre- 
quencies u! are also discrete and 

27rm . 
u; = ^ with TO G Z 

Integrating over the time interval < r < T and over 
space, one gets one has 

{Sp')^2k^\ak,^\^T 



{Sp'^)=2k^\ak,^fT 



{6f)^2Lo^\ak,u,\^T 



(SjSp) = ~2ku!\ak,u>fT 



{SpSp') = {SjSp') = 

Therefore the superposition of all the fluctuations dTOl l leads 
to 



Pro(jo, {ak.uj}) ^ exp 



Jo * 



4E 



ak,u 



2(7 L 

fjS + (J 2(72 



where some terms independent of jo have been forgotten (they 
will be fixed later by normalization). After integrating over 
the Gaussian fluctuations and if one replaces the sum over w 
by an integral one gets 



Pro(jo) ^ exp 



Jo t 



2ttL^ 

l<k<k„ 



E 



dujln 



2a L 

jma + jpa'kf D^k^ j^a"P 

(71) 
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son for these cut-offs is that the macroscopic fluctuation the- 
ory i65[ is valid only on hydrodynamic space and time scales. 
For X — 0{L^^) or r = 0{L^^) it has no validity at 
all, meaning that the cut-offs should satisfy /cmax < L and 

'J-'max < L . 

For large L, i.e. for large fcmax and iximax^ one can see by 
integrating over oj that only the constant term and the term 
proportional to Jq depend on the cut-offs so that 




a'' a 

•2 



5 ( fcmax ) '^inax)jo 



4 _ Jo" 



2(72 



2a ^ ' ADcr 



(72) 



Then one has 



K - {K) ~ 2L^ f [ dx[{6p' 
Jo Jo 



)-5p{x,Tf 



Then one can proceed as above ( I69II74I ) and get, up to terms 
constant or proportional to s, in the exponential 



^-s{K-(K))x 



ill 
2a L 



L^-^ V o"^ a 2a^ 



The rest of the calculation is the same as (1721174b . with a max- 
imum over jo achieved at jo = 0, and one finally gets 



(76) 



which is exactly (fT4l i. 



where we have used the definition (fTSI l of T. 

If the averaged rescaled current is jo over a macroscopic 
time T, the sum of the microscopic flux over all the bonds is 



Q — TL^jo = tjo. Thus as limt 



-a (see (EB) 



one can determine the cut-off dependent constants and get 



Pro(jo) exp 



2aL^ L L2 



IQD^a 



(73) 

where J- is defined in (fTSI l. This becomes, at order l/L"^, 
using the fact that iPq{s) — maxj,, [— jOs + t^^ InPro(jo)] 



HQ' 



2t 



16L»2 



(74) 



This formula is in principle valid for arbitrary diffusive sys- 
tems, i.e. for arbitrary functions a{p) and D{p). As a — 
2p{l - p), D = 1, a" ^ -4 for the SSEP this leads to the 
announced result ( |23|16t . 

For a general diffusive system the expressions of the cumu- 
lants ( l22b would therefore become 



lim 

t — *oo 



t 



— B2n-2 



(2n) 



n\ (n — 1)! 



D 



8D^ 



(75) 

where a{p) and D{p) are the two functions which appear in 
I and the -B„'s are the Bernoulli numbers. 



C. Calculation of ipg in tlie case of a weak asymmetry 

One can also repeat the above calculation in the case of 
weakly driven systems, i.e. for systems where there is an ad- 
ditional driving force of strength 1/L. This would in partic- 
ular be the case for the weakly asymmetric exclusion process 
(WASEP) lfT6ll for which the hopping rates to the right and to 
the left are respectively exp and exp(— ^). 

For such systems, ( |65] ) becomes 

Fro{{p{x,T),j{x,T)}) 
Jo Jo 2CT(p(a;,r)) 

(77) 



exp 



Following exactly the same steps as before, one gets an ad- 
ditional term ^ Sp^ in ( [69] l, everything else remaining the 
same. Then (l73T l becomes in this case: 



Pro(jo) ~ exp 



DT 



{jo-vaf{L-l) t 
2aL^ L 



IQD^a 



(78) 



where we have adjusted as in ( |73] ) the terms linear and 
quadratic in jo which are cut-off dependent. 



B. Calculation of ipK for the SSEP and derivation of fl4l ) 
To obtain ( fT4b . one can first write the activity K as 

K^2L^( dr ( dxp{x,T){l- p{x,t)). 
Jo Jo 



D. Phase transitions 



The function J^{u) becomes singular as u ^ ^ (see ( fT6b ). 
For systems for which a" < 0, this implies the occurrence 
of a phase transition in the expression ( |76l ) of V'/f (s) in or in 
the large deviation function jTSl ) of the current in the case of a 
weak asymmetry. These phase transitions are exactly the same 
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as the one discussed in fgl [lOl [161 [17|] : beyond the transition 
the system does not fluctuate anymore about a flat density pro- 
file, but the profile becomes deformed on a macroscopic scale. 

For systems such as the Kipnis Marchioro Presutti model 
lISTi [58J which have a" > 0, a similar phase transition occurs 
in ij}Q even in absence of a weak asymmetry. 



V. CONCLUSION 

In the present paper we have obtained exact expressions 
(II2I2II 1 of the first cumulants of the activity K and of the in- 
tegrated current Q for the SSEP. In the large L limit, these 
cumulants take scaling forms ( I13I22| |. 

We have shown in section|lIl]that these scaling forms can be 
understood starting from the Bethe ansatz equations ( 127159b . 
by calculating the leading finite size corrections. These finite 
size corrections are similar to the ones calculated recently for 
spin chains in the context of quantum strings [29.,, 30.1 . 

We have also shown in section |IV] that they can also be 
understood starting from the macroscopic fluctuation theory 
(l65T l of Bertini, De Sole, Gabrielli, Jona-Lasinio and Landim. 
This enabled us to extend ( I74I73I75T | our results for the SSEP 
to arbitrary diffusive systems and to see that the occurrence 
of phase transitions can be predicted from the scaling form 
of the cumulants of the current. In order to better understand 
these phase transitions it might be interesting to characterize 
the eigenstate of the s-dependent evolution operator by, e.g., 
determining correlation functions in those states. 

We have discussed here systems governed by diffusive dy- 
namics with a single conserved field. How the universal scal- 
ing forms would be modified for systems with several con- 
served fields is an interesting open question. 

We thank N. Gromov, H.J. Hilhorst, V. Kazakov, K. 
Mallick, S. Prohlac, H. Spohn, P Vieira, R.K.P Zia, for 
several useful discussions. This work was supported by 
the French Ministi-y of Education through an ANR-05-JCJC- 
44482 grant and LHMSHE. 

APPENDIX A: SEVERAL REPRESENTATIONS OF THE 
FUNCTION T 

In this appendix we show the equivalence between several rep- 
resentations ( fT6l 117152b of the function T defined in (fTSI l 



T {u) = -4 ^ [n7rv/rj2^r2~~2 

n>l 

To do so consider the integral I 



(79) 



I— / dy coth(u\/ 1 — y"^) 

I" J-i 



Then by using the fact that 

coth z = - 



2z 



n=l 



and by integrating over y, one gets 

,3 /•! 



/= / y"^ dy coi\i{u\/ 1 - y^) 

Ti" J-i 



(80) 



ri>l 



u 

"2" 



This establishes ( |52] |. Now as 



-X" = 1- 



X 



e-^ — i ^ — ' 71! 

n>0 



2 12 720 30240 



(81) 



which is simply the definition of the Bernoulli numbers 
(so that B2 = ^, -84 = ^■JH' ^6 = i^: ■••)^ one can show that 



2fe-2^2fc-3 ^1k-2 



Therefore 



fe>2 



2y3 /•! 

/ y2dycoth(uVl - y^) 



(2fc- 2) 



(82) 



2u 



2 /.I 



■E 

fc>2 



TT (2/C-2)! 



y^[\-y^) ^ ' dy 



k>2 



T{k)T{k + l) 



(83) 



Comparing dSOl l and ( |83T l. one gets 



k>2 



T{k)T{k + l) 



{-2uf 



3 45 378 2700 
(84) 

so that ( [TtI i and ( fTsT i are consistent with ( fTSb . 

For large negative u, one gets, by replacing in ( |79] | the sum 
over n by an integral. 



T{u) 



27/2(_^)3/2 

37r 



(85) 



APPENDIX B: CALCULATION OF THE TWO SUMS 
APPEARING IN ^ 

In this appendix we calculate the two sums which appear in 
( |35] l when 5^0 and L ^ 00 keeping LS fixed. 

The first sum in ( [35^ : 

If the ki are distributed according to a density g{k) on the 
real axis one can write that 



L I g{k')dk' = n 



(86) 
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Therefore for n fixed and large L, one has 
so that 



29' {h) 



n v?g'{ki) 



g{k,)L 2g(fc,)3L2 



(87) 



Replacing kj by expression dSTl i into the first sum in ( |35] ) one 
gets 

2ff-(fc,)(l-fcf) 



+ (1 - klYgihfLH-^ 



Using the fact that for no > 1 (and b < 0(1)) 



"0 ^ 
E 1 



1 TT 

+ — coth7r& (88) 



n2 + 62 262 26 
the first sum in ( l35] t can be replaced by 



2ki 



9'ih)ii-kf) 



- 1 



9{h) 

7r(l - ki)g{ki)L5cot\i['K(l - kf)g{h)LS] 

(89) 



The second sum in i 

Let us consider the following integral. 



I^V g{k')Ak' 



ki k 

We are now going to compare this integral with the sum 

1 k^kj 



(90) 



kj k ■j 



jif[i-no,i+no] 

We assume ( l86l l that the kj are given by 

L I 9{q)dq = j -a 
J-0 

and for the moment a is arbitrary. Therefore 

1 



kj+i - kj 



g{h)L 



(91) 



(92) 



One can decompose the integral / as 



/■ '+"0 1 _ t o 

/ = P / gmq\-^ 



E / 5('7)dg^3^ + E / s{q)A<iT 



1 - fcjQ 



3 = 1 ""^ 



9iq)dq 



q 

h~q Jk„ ki~q 



j=i+no '^1 

1 - hq 



(93) 



As kj+i — kj is small and of order 1/L and because of (19 11921 ) 



kj q L ki kj 

g{kj){kj+i — kjY d /I — kik^ 



1 1 k'i k j 



dkj \^ kj^ kj 
1 d / 1 — kikj 



L ki — kj 2L^g{kj) dkj \ ki ~ k 
1 1- fcifcj+i 1 

L ki — kj j^i 2L^g{kjjf-i) dfcj+i \ ki — k 



d ( \ — kikj+i 



(94) 



Therefore using ( |94b in the sum l<j<i — fiQ — 1 and 
(|94] | in the sum i + no<j<N — 1, one can rewrite (|93] | as 



fc,+,>o 1 _ /r-n 1 '-"V^ ^ -k k- 



^ 1 — kikj 1 * "'^ 1 d f 1 — kikj 

0T2 'jr. 



- y 

L . ^-^ . fc,: — fcj 2L2 ^fcj \ ki ~ k 



i=i+no + l 



\ ^ I d 

9T2 TwTTi' nr. 



1 d f \. k'i k j 



I 9{^j) dkj \ ki kj 



2L2 

j=i+no + 

9{qm-r + / 9{q)Aq-r 

-e h -q JkM f^i-i 



This becomes 



'+"0 1 - fc,(Z 1 * ^ 1 1 - k,.k 

9iq)dq- — r + 7 2^ 



kj k -j 



1 ^ 



j=i+na + 



^ kj — kj 



2L 



1 kj k- 



'if^i—no — l 



7— no — 1 



1 - kikj^np+i _ 1 - hki ^ 1 - kjkN 
ki ki-\-rLQ-\-i ki ki ki k]\[ 



9{q)M 



1 - hq 

ki - q 



9{q)dq 



1 - kiq 
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which can be rewritten as 



ki-r. 



4 E 



TV 



1 k j 



L ^-^ ki 

j=i+no + l 



1 

2L 



1 kj k 



1 kiki- 



"0-1 



+ 



1 — kiki 
k\ 

1 - kikN 



1 

1 

^2L 



ki 



i-z— riQ — 1 
I 

5(9)dg 



(95) 



From (|87] ) one can show that 

l-hq 2hnQ kiUoil - kf)g'{ki) 



ki. 



Lg{ki) 



(96) 



and that 



i—no-l 



(97) 

Lastly because one expects the symmetry kj = — /cAr+i-j and 
because L J^g g{q)dq ~ N, one gets that a = 1/2 in ( |9T| l and 
therefore the last two terms of ( |95l l vanish. 
Then using ( |96l97t into ( |95] l. one gets that 



— no — 1 



1 k^kj 

L kj k .) 

i=i * 



7 E 



- y 



1 k-i k j 



(98) 



/ - - I 2fc,: 



gih) 



no 



where the integral / is defined in ( l90l ). Lastly using the fact 
that g{k) — g{~k), one can rewrite the integral / in ( |90l ) as 



■P 



g{k')Ak' 



l-k" 



(99) 



so that ( I98I I becomes 

i-rio-l 



T 1^ T~k~ ^ L . k,-k 



1 /l^ ^ ^ 



2fc. - (1 - kj) 



j=i+na + l 



g(k')Ak' 



1 



U'2 



2^9'{ki) 



no 



1 
2 



(100) 



Note that ( |9T] l is not accurate for i close to 1 or N, i.e. near 
the singularities of g{k). A more detailed analysis of these 
two neighborhoods would only contribute to higher orders in 
the 1/L expansion 1,30.1 . 



APPENDIX C: SOLUTION OF THE AIRFOIL EQUATION 

In this appendix we show, in the spirit of |59], that the solution 

(j){x) of 



f{x)=V j dy 



\ <f>iy) 



C 



1 



dy 



1 - f{y) 



I _l y 1 — x^ y — X 

This solution is used to obtain ( |4TI ) as the solution of 
Let us choose 



(101) 



(102) 



(103) 



Then for x ^ [^1,1] and a ^ [^1; 1] one can see using (111 II ) 

»i 



dy^ 

1 y-x 



and therefore 

m^vf dyA 

7-1 y-x 



1 



- 1 



(104) 



(105) 



Now the following integral of this function f{x) can be com- 
puted (using ( 11091111b for for x ^ [-1, 1] 



dy 

1 y~x 



(106) 



a — X 



1 + \/x^ 



1 - X 



so that 



9 



\y^^^m 

1 y-x 



a — X 



= a — ^/a^ — 1 — 



- 1 - x 

1 - x2 

a — a; 
(107) 



Comparing with ( I103l l we see that 
1 



dj/ 



1 - f{y) _ a - Va'^ - I \/l - x"^ 
1 — x"^ y — X 



VT-- 

a — vc 



(108) 



Therefore ( I102l i is the solution of ( IIOII ) with a constant G 
which depends through a on (t>{x) when one chooses ( 11031 ) 
for 0(x). 

As the inversion formula ( 11021 ) is valid for arbitrary a, it 
would also be valid when f{x) is any polynomial in x, and 
as the polynomials are dense in the set of continuous func- 
tions on (—1, 1), one can consider that ( |101|102| ) are valid for 
"arbitrary functions" f{x). 
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APPENDIX D: USEFUL INTEGRALS 



and therefore for any function F{y) 



In this appendix we list a few integrals which are used in var- 
ious places of the paper. 

First for x ^ [—1,1] one has 



1 



dy ^ \/x^ -1-x (109) 



so that 



TT J -I y-x 



dy ~ —X 



(110) 



As a consequence of (II 1 11 1 one has for x ^ [—1,1] and 
a 4 1-1,1] 



(111) 



TT J_i y — X y — a a — X a — x 
and thus for x E [—1,1] and a ^ [—1,1] 



vr 7_2^ y^x y — a a — x 



12) 



One can also show that 



1 

and that for y ^ [—1, 1] 

da; 1 



(113) 



-1 vi^l? a; vyrrr 

As a consequence of ( II 1211 14l i. one has 

r^^r^F(,).o (115) 

J-iVl~x^ J-iV-x 
for an arbitrary function F{y) as it is valid for any polynomial 

For 9 < 1 one can show using dl 14b that 

(116) 



dx 



Li {i-e^x^)VT^~^ 

one can also show 



' VT-x 



_1 (1- 6*2x2) 

and that 



dx ~ n — (117) 



02 



.1 l-6'2y2 

and for y ^ [—1, 1 
1 



d'y = T^\-^-7^ 



1 1 VT~e^ 



6*4 26*2 6*4 



dy 



(118) 



1 (l-6l2x2)yr^^ y-x (1 _ 612^2)^^^2—1 



776*2 y 



(1 - 612y2)yr3^ 

(119) 



dx 



dy 



-1 (1 - 6'2x2)Vl^:^ J-iy-a: 

2 yF(y) 



dy 



(120) 



APPENDIX E: BETHE ANSATZ CALCULATION FOR 
THE CURRENT LARGE DEVIATION FUNCTION i/'q(s) 

This appendix describes how a Bethe ansatz calculation of 
ipqis) similar to the one conducted for ipK can be imple- 
mented. The operator Wg whose largest eigenvalue is tpq 
reads, in the spin language aheady used in ( 1251) . 



L r 



- e V+(7,+i + eV, (7+1 



(121) 

The Bethe ansatz equation analogous to dZTl i take the form 



N r 



n 



1 - 2e-^C. 



-2s 



CO 



l-2e--0+e-2-C<0 



(122) 



In terms of the Q 's, we have that 



(114) V'q(s) =-2A^ + e ^[Ci + ... + Civ] 



1 



1 

_ (123) 
Kim 152] has studied the spectrum of 7^ = -WQ/(coshs/2) 
by means of a Bethe ansatz calculation: in the notations of his 
equation (1), the parameters A and S are given by 



A = 



1 



coshs ' 



S = tanh s 



(124) 



but unfortunately his results do not apply to our particu- 
lar case, which turns out to correspond to a critical point 
of the related six-vertex model. The defining parameters of 
the latter, denoted by A, H and ly, are related to Kim's by 
A = A/cosh(2iJ), S = tanh(2iJ), A coshiy. Thus, in 
terms of our original parameters, we get that 



A = 1 



2H 







(125) 



a limiting case explicitly excluded by Kim which lies at the 
critical point of the six-vertex model. 

We choose to write that Q — ^-^i^if^i+'^^p) ^ The two 
main differences with the calculation of ijjK is that the C^'s 
dependence in s is different. We have also shifted them by 2ip 
for convenience. Just as was the case previously, the fc^'s will 
be densely distributed on a connected curve C of the complex 
plane that is invariant upon complex conjugation. Given 
that the equations for the Cj's are invariant under complex 
conjugation, we expect the contour C to be symmetric with 
respect to the vertical axis in the complex k plane. We shall 
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denote the end points of C by —0* and 6. 

Given that (1122b becomes 
1 ^ 

—i{ki+i2p) = — U{ki,kj), where 



The four unknowns C, 0, r and /i are determined by requiring 
that 00 remains finite as a; — > 1 and as a; — > —r, and by noting 
that by definition 



(126) 



U (ki, kj) = — In 



1 - 2e-^C. 



'GO 



1 - 2e-'^0 



-2s 



00 



(132) 



for |i — j| S> 1 we expect that 



2i(fci + ia){kj + za) 

A^2 ~ ^ 7 



while for i — j of order 1 , s will be over order 1/L and ki — kj 
as well. We define g{k) as the root density along contour C, 
so that 



while (/)o must verify the self-consistency equation 
h = i{a + l)/2 + J dy4ia{y){y + lajBY^. After explicitly 
evaluating the latter integral and that appearing in (1132b . we 
, a = 2p - 1 (127) arrive at r = 1, /i = and 47rC = B = 2^p(l - p), which 

leads to 4><i{x) = — . Up to a sign, this is exactly the 
same function as that found in the study of K, and this is the 
same end point 9 — 2^ p(\ — p) for the contour on which 
the fcj's lie. 



L 



g{k)dk ^j-i 



We may now simplify ( |129b into 



(note that g{k) is in general complex but along the contour 
g{k)dk is real). If kj and ki are n roots apart, we have that 

^0 ^^i = JiJrjL ~ 2g(fc03L2 + • ■ •■ Expanding U at fixed 
sL in powers of leads to 



a + 1 



1 n~ig{h){h+ia)^sL 
s n + ig[ki){ki + lay sL 



1 a;(a; + ia/OY' 



J y-x 

{\6^/l-x'^{sL) 12] coth[0v/l - a;2(.sL) /2] 

(133) 



+ la) 2 + 



g'{ki){ki + ig) 



n2 + + ia)2(sL)] 

(128) 

Equations ( 1127b and ( 1128b play a role analogous to ( l33T l and 
( |32] | in the study of After using the methods of appendices 
B and C we arrive at the following equation for g which we 
express in terms of 0(x) — {6x + ia)'^g{dx) and r — 9* /d: 



whose solution reads (p{x) = (poix) + S<j){x), 



S(l){x) 



5C 



2Sha 



1 



a;2 27r\/l — 



1 



dy 



(134) 



a + 1 



^ ( 



2^ .^;j^ 0(y)-(2/-^)(y + »t)^V(y) 



aN^2 cf)'{x) 



'K(f){x){sL)\ cot\i['K(j){x){sL)\ 
(129) 



7r2 a/1 - a;2 ' 7 " 2/ - a; 
We have denoted by 5F{x) the function 

e x{x + ia/OY 



5F{y) 



5F{x) 



2L 



1 — X 



Let us denote (j>Q (x) the solution of the above equation, in the 
L ^ CO limit 

ex/2 + h^V [ dy^^ (130) 

J^r y~x 

where h = i{a + l)/2 + J dy(po{y){y + is a density- 

dependent constant to be determined. The general solution of 
( fnob can be written (see (11011102b ) as 



X (^[9^/1 - x^{sL)/2] coth[6'v/l - a;2(sL)/2]) 



x{x + ia/9y 



2L 



l-a;2 



y^(^dsLY{l-xy/^ + l 
P- 



(135) 



The new constants SC and Sh are determined by 



(l)o{x) 



C 



(r + 1) 



/ 



and Sh ^ f -, — ; 



After performing 



^(1 - x){r + x) 167r^(l - x)ir + x) 
9x^ 2h{r -1 + x) 



explicit integrations along the lines of appendix D, we obtain 
the final result through the following equality 



27rv/(l - x){r + x) 27rv'(l - x){r + x) 

9x(r - 1) 
47rY/(l - x){r + x) 



(131) 



da;- 



1 



da;0(a;) 



-s2 + s2jc'6/7r 



(136) 



y-x 
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where 



32 /.I 



27rL 



da; X 



p>2 ^' 

1 



(137) 



After noting that, as before, we have 



it only remains to substitute the value of SC into ( I136l l. This 
allows us to conclude that 



V;q(s) = -js'iL + 1) + L-^Ti-L's'e'/S) (139) 



y-x 



SF{y) = (138) which is the announced result of dMT l. 
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